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Mean first-passage time for random walks in general graphs with a deep trap 
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We provide an explicit formula for the global mean first-passage time (GMFPT) for random walks in a general 
graph with a perfect trap fixed at an arbitrary node, where GMFPT is the average of mean first-passage time 
to the trap over all starting nodes in the whole graph. The formula is expressed in terms of eigenvalues and 
eigenvectors of Laplacian matrix for the graph. We then use the formula to deduce a tight lower bound for 
the GMFPT in terms of only the numbers of nodes and edges, as well as the degree of the trap, which can be 
achieved in both complete graphs and star graphs. We show that for a large sparse graph the leading scaling 
for this lower bound is proportional to the system size and the reciprocal of the degree for the trap node. 
Particularly, we demonstrate that for a scale-free graph of size N with a degree distribution P(d) ~ d 1 
characterized by 7, when the trap is placed on a most connected node, the dominating scaling of the lower 
bound becomes TV 1-1 / 7 , which can be reached in some scale-free graphs. Finally, we prove that the leading 
behavior of upper bounds for GMFPT on any graph is at most TV 3 that can be reached in the bar-bell 
graphs. This work provides a comprehensive understanding of previous results about trapping in various 
special graphs with a trap located at a specific location. 
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I. INTRODUCTION 

Trapping problem first introduced in the seminal work 
of Montroll in 1 is a kind of random walk in which a per- 
fect trap is located at a given position, absorbing all par- 
ticles (walkers) visiting it. Trapping process is closely 
related to various other dynamical processes in diverse 
complex systems, including lighting harvesting in an- 
tenna systems^ 3 -, energy or exciton transport in poly- 
mer systems^, charge transport in amorphous solids 5 , 
phtoton-harvesting in cells 6 , and so forth. In view of 
its practical significance, it is of great interest to study 
trapping problem in a large variety of complex systems. 

The highly desirable quantity of trapping problem is 
mean first-passage time (MFPT)2. The MFPT for a node 
i to the trap is the expected time taken by a particle 
starting from i to first arrive at the trap. The aver- 
age of MFPT over all source nodes in the system other 
than the trap is often called global mean first-passage 
time (GMFPT), which is very helpful for understand- 
ing the kinetics of the whole trapping process and pro- 
vides a useful indicator for the trapping efficiency. It 
is thus very important to evaluate GMFPT on different 
systems. Thus far, there has been considerable inter- 
est in computing GMFPT on various systems in order 
to uncover the effect of system structure on the behavior 
of GMFPT, such as regular lattices with different dimen- 
sionsi*2i£, the Sierpinski fractals^—, the T— fractal^— , 
Cayley trees^ and Vicsek fractals^ as models of polymer 
networks 1 - - — , as well as many scale-free graphs^—. 
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Although previous works have uncovered some nontriv- 
ial effects of different structural properties (e.g., scale- free 
behavior. 3 !, modularity 3 ^, and fractality 3 ^) on GMFPT, 
they focused on special graphs with the trap fixed on 
a special node. So far, the determination of GMFPT on 
general graphs with a trap at an arbitrary position is still 
open, and no quantitative formula has been obtained for 
GMFPT on a general graph. On the other hand, different 
topological features are often dependent on one another, 
e.g., in real systems fractal characteristic is frequently ac- 
companied by a modular structure^. It is thus difficult 
to discern the roles of different structural properties on 
GMFPT. Sometimes, one may confuse and even distort 
the roles of different topologies. 

In this paper, we perform an in-depth research about 
trapping problem on a general graph with the deep trap 
placed at an arbitrary node. We derive an implicit ex- 
pression for GMFPT, in terms of eigenvalues and their 
corresponding eigenvectors of Lapacian matrix associated 
with the graph. On the basis of the obtained formula, we 
then provide a lower bound of GMFPT, and show that 
it can be achieved in two graphs, complete graphs and 
star graphs. In addition, for large sparse systems we find 
that the leading term of the lower bound is proportional 
to the system size and reciprocal of the degree of the trap 
node, irrespective of any individual structure. Particu- 
larly, we show that in a scale-free graph with N nodes 
obeying degree distribution P(d) ~ d -7 , where the ex- 
ponent 7 characterizes the extent of inhomogeneity, the 
dominating scaling of the lower bound of GMFPT to a 
most-connected node is A^ 1-1 / 7 . Also, by making a com- 
prehensive analysis of previous studies, we show that our 
work can unify existing results about GMFPT obtained 
for various scale-free graphs. Finally, we show that the 
GMFPT of any graph is at most N 3 . 
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II. GLOBAL MEAN FIRST-PASSAGE TIME ON 
GENERAL GRAPHS 

In this section we study trapping problem on a general 
graph G with TV nodes and E edges, which is a particular 
random walk performed on the graph with a single trap 
located at an arbitrary node. The random walk consid- 
ered here is a simple discrete-time random wal k 41 i 42 . At 
each time step, the walker jumps from its current posi- 
tion to any of its neighbors with the same probability. 
Let Tij denote the MFPT from node i to node j. For 
convenience, let j be the trap node, and let Tj denote 
the GMFPT to node j. Then, by definition Tj, the most 
important quantity for the trapping problem, is given by 



Tj 



N 



(1) 



In the sequel, we will use the connection 4 ^ between re- 
sistance distance and MFPT to derive an explicit formula 
for Tj , based on which we will provide a lower bound for 
Tj, as well as the leading scaling for this lower bound 
when the graph is sparse and very large. Moreover, we 
will give an upper bound for Tj. 

A. Explicit expression for global mean first-passage time 

Let G represent the corresponding resistor network 4 ^ 
of graph G, which is obtain from G by arranging a unit 
resistor along every edge of G. Let R a b be the effective 
resistance, i.e., resistance distance 45 , between any pair of 
nodes a and b in the resistor network, which is defined 
as the voltage when a unit current enters one node and 
leaves the other. Then, based on the established relation 
governing MFPT and effective resistance^, we have 



T, , 



Riz) 



(2) 



where d z is the degree of node z in graph G. 

Thus, in order to determine Tj, one can alternatively 
evaluate the terms about effective resistance on the right- 
hand side of Eq. @. It is known that effective resistance 
between any two nodes can be determined by eigenvalues 
and eigenvectors of Laplacian spectrum of the original 
graph. For graph G with size N , its Laplacian matrix 
L is an N x N matrix, whose elements are defined 
as follows: the non-diagonal entry Zy = —I, if nodes i 
and j are connected by a link, otherwise hj = 0; while 
the diagonal entry la = d{. Let Ai, A2, A3, \jq 
be the N eigenvalues of L, rearranged in an increasing 
order as = Ai < A2 < • ■ ■ < Ajv, and let /j,x, (1%, ■ ■ -/xjv 
be the corresponding mutually orthogonal eigenvectors 
of unit length, where \ii = (nn, fia, ■ ■ ■ Hin) J ■ Then the 
effective resistance Rij is given by2& 



Ri 



1 



N 

k=2 K 



HkjY 



(3) 



Inserting Eqs. ([2]) and ([3]) into Eq. ([TJ), we have 
I 1 N 1 



k=2 



(4) 



Equation (|4]) can be simplified by using the properties 
for related quantities of matrix L listed below. According 
to the spectral graph theory 4 ^, we have: 

(i) the eigenvectors of matrix L satisfy 



li]Hi = Sij , 



(5) 



where Sij is the Kronecker delta function defined as: Sij — 
1 if i is equal to j, and 5y=0 otherwise. 

(ii) the entries of eigenvector /ii corresponding to Ai = 
are all equal to each other, namely = \A/V /N holds 
for 1 < i < N. Then, for 1< k < N, 



^Hki = 0. 



(6) 



(iii) the entry lij of matrix L has the following spectral 
representation: 



N 



hj — y ' ^k^ki^kj 
fc=l 



which means 



di = A fc 



2 



(7) 



(8) 



k=l 



The above introduced properties of matrix L are very 
useful for the following derivation. We first use them to 
simplify Eq. ([J} that can be recast as 



(Tj) = w^r E ^ E E -EE d *^kj 

k—2 i Z i z 



(9) 



From the above intermediary results, the four terms in 
the square brackets of Eq. (|9]) can be sequentially calcu- 
lated as follows: 



Y,J2 d ^l 3 =Nx2Ex,i 2 kj 



(10) 



and 



\] }^ dziJ, k iLi k j = 2E x /Zfcj jiki = 0; (11) 

i z i 

E E d ^^kj^kz ~ N x fi kj d z [J.kz , (12) 



^^dzfikiUkz = ^2^ki^2d z ^ kz = 0, (13) 

i z i z 
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where the facts E — ~}2, z d z /2 and A = 1 are used- 
Plugging Eqs. (I10H13[) into Eq. ©, we arrive at the ex- 
plicit expression for GMFPT on a general graph with the 
immobile trap fixed at an arbitrary node j, given by 



which is a main result of this work. It is not obvious that 
Eq. (|14p depends on any individual structural feature. 



B. Lower bound for global mean first-passage time 



T, 



A 



A 



fc=2 fe V z J | 



After obtaining the exact formula for GMFPT, we next 
apply it to derive a lower bound for Tj. By Cauchy's 
inequality, one has 



N I ( 

A=2 fe V z 



Afe IE x filj - (J-kj 



E d 



z fJ k z 









> 







' N ( 

,fc=2 \ 



E d 



(15) 



Considering Eqs. (0 and dHJ and Ai = 0, we have 
Y2 x k ( x /4j ~ Mfcj E d ^ fcz ) 

fc=2 V z J 

= 2E ^2 ^k^kj ~ E ^ z E ^k^kj^kz 
k z k 

= 2Edj -Y,dzl 3 z- (16) 

In addition, 
N ( 

k = 2 \ z 

= 2E^2 Aj - - E dz E M fe jMfcz + E dz/il2 



Equation ([5]) shows that 



A' 



E ^kjfJ'kz — 



k=l 

which leads to 



1, if j = z, 
0, otherwise, 



}^d z }^jj lk jfikz = d 3 



z k 



and 



Mi? E dz ^ lz 



2E 
~N~' 



(17) 



(18) 



(19) 



(20) 



Therefore, Eq. (fTTj) can be simplified to 

^kj — Mfcj 
fc=2 \ z j 



2E - 2E x — - dj + 2E x — = 2E - cL. (21) 
A J AT 



Combining the above obtained expressions, we have 
J ~ N — 1 2E x dj — J2 Z dzljz ' 



(22) 



where Zj 2 is the entry of Laplacian matrix L as defined 



J-(2£^L - ^kjJ2z d ^^kz) and 6 fc = 



before. 
Let a/e = 

\ k (2En 2 kj - HkjJ^z^^kz)- In Eq. fl22J|, the equality 
holds if and only if 02/62 = 03/63 = • • • = clm /&jv- We 
distinguish two cases. If 2E^ — fikj J2 Z dzHkz 7^ for 
all 2 < A; < A, then the equality of Eq. (g2J holds if 
and only if A2 = A3 = • • • = Xn- In this case, only 
when G is isomorphic to the complete graph, the equal- 
ity holds 48 . For the other case that there exists some 
k such that 2E/i1j — Hkj J2 Z dz^kz = 0, then the equal- 
ity holds if for those k with 22?/iL — Hkj J2 Z d z /ikz 7^ 0, 
the corresponding eigenvalues Afc are equal to each other. 
For example, the equality holds for the star graph, since 
for the star graph the term 2E[}^.^ — [ikj J2z dz^kz in 
Eq. (|15[) is non-zero only when k — A, while for other k 
(k = 2, 3, • • • , N - 1), 2Enlj - fikj Y.z dzVkz is exactly 
zero. 

It is easy to prove that the term ^ z d z lj z reaches its 
minimal value when node j is linked to all other nodes 
in the graph excluding node j itself, that is, 



dj 



jz> d < 



(2E-d j ) 



(23) 



Therefore, 
N 



(2E-d J 

3 ~ A - 1 2E x d 3 ; - d? + 



A 2E- d-i 



2E - d., 



A — 1 dj + 1 

(24) 

In this way, we have obtained a lower bound for GMFPT 
in terms of the number of nodes, the number of edges, 
and the degree of trap node, but independent of other 
structural parameters. 

The lower bound for GMFPT given in Eq. (p4|) is sharp 
since it can be achieved in some graphs. For example, it 
has been reported that for the complete graph with A 
nodes, the MFPT between any node pair i and j is ex- 
actly A — 1—. Thus, for trapping on the complete graph 
with the single trap positioned at an arbitrary node, the 
GMFPT is A — 1, which is equal to the lower bound pro- 
vided by Eq. (|24l) . Again for instance, it can be proved 
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that for trapping on the star graphs with the trap being 
placed on the center node, Eq. (IM1) can also be reached. 

The obtained lower bound for GMFPT, i.e., Eq. 
provides important and useful information of a graph 
since it establishes a range of GMFPT in terms of the 
simple graph parameters. Moreover, although the lower 
bound can only be achieved in some very few cases, as we 
will see that its leading scaling can reconcile previous re- 
sults on GMFPT for random walks on different networks, 
especially on scale-free networks. 

Extensive empirical works have shown that most real 
networked systems are sparse in the sense that their aver- 
age node degree (d) = is a small constant 5 -^— . Equa- 
tion (|24l) shows that in sparse graphs, when trap node is 
placed on a node with degree d, the lower bound becomes 



T d > 



N N(d) - d 
N - 1 d+1 



(25) 



C. Scaling of upper bounds for global mean first-passage 
time 

We proceed to provide an upper bound for the GMFPT 
on a general graph. For this case, it is difficult to provide 
an exact expression as that of the lower bound, but we 
can give a scaling of upper bound. To this end, we recall 
another quantity, i.e., cover time, for random walks. The 
cover time for a node i denoted by Ci is the expected 
time needed for the walker starting from i to first visit 
all the nodes in the graph. It has been established— that 
for an arbitrary graph, the upper bound of cover time Cj 
grows with the system size N as iV 3 . Note that in any 
graph, it is evident that Ty < Ci, implying that N 3 is 
also an upper bound for GMFPT of the graph. It is easy 
to proved that this upper bound can be reached in the 
bar-bell graph of N nodes, which consists of two cliques 
each of size AT/3, connected by a path of length N/3. 



In the limit of large network size N, the leading term of 
the lower bound is 

T d >^, (26) 

which is proportional to the system N and the inverse 
degree of the trap node. We note that this dominat- 
ing term has been previously reported in 3 ^, deduced by 
another technique, but the definition of GMFPT given 
ir* 3 ^ is different from the one here that is more frequently 
adopted in the literature. 

In addition to the sparseness, most real-life networks 
also display scale-free properties 3 - 7 - with their degree dis- 
tributions P(d) obeying a power-law form P(d) ~ d~ 7 
characterized by the exponent 7, which suggests that 
there exist large-degree nodes (hubs) in these networks 
whose degree d max satisfying d max ~ A rl ^ 7 ^ 1 ^ 53 ' 54 . 
These nodes have a strong effect on dynamical processes 
occurring on graphs 55 , including the trapping process dis- 
cussed here. For example, when a trap is located at a hub 
in the pseudofractal web^ 7 ., the GMFPT exhibits a sub- 
linear scaling with system size N as N^ . Such sublinear 
scaling underlies a high trapping efficiency which was ev- 
idenced also i n 28 i 30 i 31 . Equation (j2"6")l implies that when 
a trap is fixed at a hub with degree <i max on a scale-free 
graph, the scaling of the lower bound for GMFPT is , 

behaving as jV 1-1 ^ 7-1 ) in large systems. This minimal 
scaling is the possible highest trapping efficiency that can 
be achieved in some hierarchica l 30 ' 33 ' 36 and modula r 31 ' 35 
scale- free graphs. This implies that, in the future, look- 
ing for special topologies exhibiting an efficiency higher 
and higher may be a kind of obsolete problem. On the 
other hand, besides the sublinear scaling, the GMFPT 
on scale-free graphs can also display linea r 56 ' 57 even su- 
perlinea r 58 ' 59 scaling with the system size, all of which 
are encompassed in the framework presented here. 



III. CONCLUSIONS 

We have studied the trapping problem on a general 
graph, which is a special random walk with a trap lo- 
cated at a certain node. Using the spectral graph theory, 
we have obtained an explicit solution to GMFPT, which 
is expressed in terms of the eigenvalues and eigenvectors 
of Laplacian matrix of the graph. Based on this result 
we have further provided a lower bound for GMFPT that 
can be achieved in complete graphs and star graphs. The 
obtained lower bound implies that in large sparse systems 
its leading term grows proportionally with the system size 
and the reciprocal of the degree of the target node. Par- 
ticularly, we have demonstrated that for trapping on a 
scale-free graph with the immobile trap being positioned 
on a hub node, the minimal scaling of the lower bound 
can be expressed by the exponent degree distribution 7 
characterizing the inhomogeneity of the graph, which can 
be reached in some special scale-free graphs. At last, we 
have deduced an upper bound for the GMFPT. Our work 
provides a broader view of previous researches for trap- 
ping on diverse graphs and sheds light on some aspects 
related to the trapping problem, e.g., lighting harvesting. 
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